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^ ■ Abstract 

In this paper the problem of strong solvability of the incompressible Navier-Stokes equations 

m . 

(INSE) is revisited, with the goal of determining the minimal assumptions for the validity of 
a local existence and uniqueness theorem for the Navier-Stokes fluid fields (solutions of INSE). 

■ Emphasis is placed on fluid fields which, together with suitable derivatives, do not necessarily 
decay at infinity and hence do not belong to Sobolev spaces. For this purpose a novel approach 
based on a so-called inverse kinetic theory, recently developed by Tessarotto and Ellero, is adopted. 
This requires the construction of a suitable kinetic equation, advancing in time a suitably smooth 

■ kinetic distribution function and providing exactly, as its moment equations, the complete set of 



fluid equations. In turn, by proper definition of the kinetic equation, this permits the introduction 



in 

of the so-called Navier-Stokes dynamical system, i.e., the dynamical system which advances in time 

00 ' 

self-consistently the Navier-Stokes fluid fields. Investigation of the properties of this dynamical 



system is crucial for the establishment of an existence and uniqueness theorem for strong solutions 
of INSE. The new theorem applies both to bounded and unbounded domains and in the presence of 
generalized boundaries, represented by surfaces, curves or even sets of isolated points. In particular, 
for unbounded domains, solutions are considered, which do not necessarily vanish at infinity. Basic 
consequences for the functional setting of classical solutions are analyzed. 

PACS numbers: 47.10.ad,47.10.Fg,47.10.A- 
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I. INTRODUCTION 

The goal of this paper is to investigate a well-known theoretical issue of fluid dynamics. 
This is concerned with the existence and smoothness of solutions of the fluid equations which 
characterize an incompressible Newtonian fluid, i.e., the so-called incompressible Navier- 
Stokes equations (INSE). In particular, we will be interested in solutions which are physically 
realizable, i.e., can be identified with physical observables. In the theory of continua these 
are necessarily described by strong solutions. Several aspects of the mathematical theory of 
these fluid equations remain unsolved. This occurs, in particular, in the case of: 

• 3D external domains; 

• forcing; 

• solutions which do not necessarily decay at infinity in configuration space. 

A related problem of mathematical research involves the search of the classical dynamical 
system - here denoted as Navier-Stokes (NS) dynamical system - which uniquely advances 
in time the related fluid fields that characterize the Newtonian fluid, i.e., the fluid velocity 
V(r,t) and the fluid pressure p(r,t). In other words, if the NS dynamical system actually 
exists, it will permit to cast the complete set of fluid equations in terms of an equivalent 
(and possibly infinite) set of ordinary differential equations which define the dynamical sys- 
tem itself. From the mathematical viewpoint the reason why such a dynamical system is so 
important is that, if a theorem of existence and uniqueness can be established for the set of 
ordinary differential equations which determines the dynamical system, than it is obvious 
that it will imply the same conclusion also for INSE. For contemporary science the determi- 
nation of such a dynamical system represents not merely an unsolved intellectual challenge, 
but a fundamental prerequisite for the proper formulation of all physical theories which 
are based on the description of these fluids. These involve, for example, the understand- 
ing of the related phase-space Lagrangian fluid descriptions and the consistent formulation 
of turbulence theory and of the related approximate statistical descriptions (i.e., obtained 
introducing appropriate stochastic models able to reproduce phenomenological data), both 
essential in fluid dynamics and in applied sciences. 
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the problem 
5|, is based on the 



Surprisingly, until recently (see Tessarotto and Ellero, 2000-2005 1 , 
has remained unsolved. Its solution, discussed in detail in Refs. 3), 4, 
construction of an inverse kinetic theory (IKT) for the incompressible NS equations, i.e., 
a suitable phase-space description which provides the complete set of fluid equations in 
terms of an appropriate inverse kinetic equation. Such an equation, in particular, can be 
constructed in such a way to define uniquely a classical dynamical system which advances 
in time the fluid fields, i.e., so that they result - by construction - solutions of the same fluid 
equations. For the extension of the theory to the treatment of incompressible thermofluids 
and quantum fluids see Refs.B,H,y,3. 

In this paper we intend to point out basic consequences which can be drawn from IKT 
as far as the problem of existence and smoothness of solutions of INSE is concerned. 



A. Strong solutions of the INSE problem 

For definiteness, we shall assume that the relevant fluids {p, V,p} , i.e., respectively the 
mass density, fluid velocity and pressure describing the fluid, are defined point-wise and 
suitably smooth in an appropriate domain Qxl. Here Q (fluid domain) is an open connected 
subset of the Euclidean space M 3 , with boundary dQ and closure f2. f2 is defined as the set 
where the mass density is a constant p = p Q > 0; moreover, I is either a time interval 
/ =]tojM (with closure /= [io,ti]), or the real axis R. We assume that the fluid fields are 
continuous in Q x J, fulfill suitable initial and boundary conditions respectively at t — t Q 
and on the boundary dfl, while in the open set Qxl they satisfy the fluid equations: 

-^ p + pV-V = 0, (1) 

p_^ V + Vp-f-pV 2 V = 0, (2) 
V • V = 0, (3) 

P = Po, (4) 

which are subject to the inequalities 

P > 0, (5) 
p > 0. (6) 



3 



Here the notation is standard. Thus, p a and p > are respectively the constant mass density 
and fluid viscosity and f (r,t) the volume force density. Moreover, = J^V + V ■ V is the 
convective derivative, Eq.(CE]) denotes the so-called continuity equation, while Eqs.([2]), ([3]) and 
(|4|) are respectively the Navier-Stokes (NS), isochoricity and incompressibility equations. 
Finally, the inequalities (j5J) and (jBJ) represent the so-called physical realizability conditions 
of the fluid, which must be prescribed in order that p and p are physical observables. 

The complete set of equations CQ)-(J4]) - subject to the inequalities (j^J) and (jUJ) - are denoted 
as incompressible Navier-Stokes equations (INSE) and its solutions NS fluid fields. In the 
following it is assumed that: 1) {p, V,p} are strong solutions of the INSE problem which 
are suitably smooth so that Eqs.( fT2l and (|T3l) are identically satisfied in the whole set Qx I; 
2) the solutions of INSE are assumed to satisfy suitable initial- (Dirichlet-) boundary value 
problem (INSE problem). These are defined as follows: 

A) fluid initial conditions: the initial conditions for the fluid fields A(r,t ) = {p, ~V,p}/ rt \ 
are defined imposing 

A(v,t ) = A (r), (7) 

where t Q G I, and A D (r) = {p Q , V (r),p (r)} satisfy respectively the condition of isochoricity 
and the Poisson equation: 

V • V G = 0, (8) 
V 2 Po = -V • [V -VV ] -V ■ f(r,g. (9) 

B) fluid boundary conditions: dQ is considered for greater generality as a moving boundary. 
In particular, for all points vwif) G dfl let us assume that their velocity 

V w {y w {t),t) = jv w {t) (10) 

is a suitably smooth function of time defined in /. Then, the boundary conditions are 
obtained, in both cases, imposing the Dirichlet boundary conditions 

lim A(r,t) = A w (r w ,t), (11) 
where A w (r w ,t) = {p a , Vw,Pw}( Fw>t ) ■ 
In particular, we remark that: 
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• by taking the divergence of the NS equation (j5J) and respectively its scalar product 
with V, it follows the Poisson equation for the fluid pressure p 

V 2 p = - Po V ■ (V • VV) + V • f (12) 

and the energy equation 

P WY + V ' ^ p ~ f " ^ 2y ] = °- ( 13 ) 

• since conventionally the domain of vacuum is defined as the subdomain (of IR 3 ) in 
which both p and p vanish identically, the inequalities (jSJ) and (jUJ) provide the physical 
requirements in order that the domain Q is non-empty; 

• the boundary of the fluid domain, dQ, may generally include so-called free boundaries 
where the fluid pressure p(r,t) locally vanishes; 

• in case Q is unbounded and dQ (or a subset of dQ) is an improper surface of M 3 , the 
fluid fields A(r,t) need not decay necessarily at infinity. This means that A(r\y,t) may 
not vanish at infinity As a consequence, in such a case it follows that V(r,t) shall not 
be required to be L 2 {Q). 



B. Motivations and historical background 

The modern status of mathematical theory of NS equations is largely due to the pioneering 
work of Leray in the years 1933-34 1(3, 11, 12] and Hopf in 1950-51 [3] who reformulated 



the NS PDE's in terms of a suitable set of integral equations and introduced the concept 
of weak solutions for NS equations. These solutions have played, since their introduction, 
a major role in the mathematical research dealing with NS equations. Indeed, they are the 
only solutions which, so far, have been proven to exist for all times and without restrictions 
on the initial data, apart the requirement of a suitable functional setting for the fluid fields 
(for a review see Galdi, 2002 [371] ) - However, the fundamental problem of a unique global 
weak solution and, conversely, the possibility that uniqueness holds only locally due to the 



appearance of a local "bifurcation" phenomenon (the so-called Leray conjecture 23j) still 
remain open issues. 
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On the other hand, the distinction between strong (or classical) and weak solutions has 
important relevance also for the modelling of real fluids in the framework of continuous 
mechanics. It should be pointed out that the existence of weak solutions, rather than the 
classical ones, for NS equations may be also potentially an indication of the failure of the 
mathematical model based on the incompressible Navier-Stokes equations. In fact, it may 
involve the violation of the subsidiary conditions indicated above, in particular since a weak 
solution is generally not defined everywhere in Q x / (Doering and Gibbon, 1997 3l|). 

An important feature which characterizes weak solutions is the way in which they are 
usually manufactured, i.e., their existence and uniqueness is established. This is obtained 
replacing the NS PDE's by an infinite set of ODE's and by constructing explicitly their 
solutions, which furnish a sequence of successive approximations, the so-called Galerkin ap- 
proximates (|3l|)- The existence and uniqueness theorem for weak solutions of NS equations 
is achieved by demonstrating that the Galerkin approximates converge in a suitably weak 
sense to the solutions of the original PDE's (which is another way of justifying the name 
given to these solutions). 



In its traditional approach 
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261 ] the treatment of the existence 



and uniqueness problem in the local sense, i.e., for a finite time interval /, requires a suitable 
functional setting which depends, in particular, on the choice of the boundary conditions. 
A necessary condition for the fluid fields to be strong solutions of INSE is that they are 
one-sided continuous on boundaries which separate the configuration space. Moreover, in 
the open domain Q, since all fluid equations [i.e., Eqs.(pQ)-(jI]) as well as ( |T2l and ( |T3i) ] must 
be defined by functions which are at least continuous, i.e., C^ ' \il x /), the fluid fields and 
the volume force density f (r,t) which satisfy INSE must evidently have the native functional 
setting (NF setting). These conditions involve both the fluid fields {p, V,p} and the volume 
force f(r,t), which are required to satisfy at least 



' {ftV lP }eC(°)(fixl), 
V(r,t)EC^(Q x /), 
p(r,t)eC^°\n x /), 
f(r,t)eC( 1 '°)(fix /), 



(14) 
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as well as their initial and boundary conditions, which imply also 



Po (r)eC^(n), (15) 



and 

' V w (r w (t),t)eC^(I), 



(16) 



Pw(r w (t),t)eC^(n). 

Indeed, in validity of ffT4|) [and f|T5|) . f|T6l) ]. by invoking equations ([2])-© and imposing in 
the whole set f2 the initial condition 

P {T,t ) = P o>0, (17) 

it is obvious that Poisson equation Eq. (fl2l) implies the incompressibility and the isochoricity 
conditions ([3]) and (JH), as well as the energy equation (fT5|) and vice versa. Hence, imposing 
Eq. (jT71) . INSE can be reduced to the equivalent set given by Eqs. (Op)- (BJ and the physical 
realizability conditions ([5]), ([6]). Here = U dVl is the closure of f2 and C^(Q x I) = 

c«(fi)ncW(i). 

The customary approach to weak solutions is, instead, based on the requirement that the 
fluid fields, together with the data, belong to appropriate Sobolev spaces Q, 2, 25], 



endowed with a suitable scalar product, usually defined by the integral on the configuration 
space Q 

(V(r,t), V a (r,*)) = / drV(r,*) • V 1 (r,t) ) (18) 
Jo, 

being V(r,t), Vi(r,t) two independent velocity fields, and with norm 

||V(r,f)||;= [ drJ2 (D a V,D a V). (19) 

Jn a=0,l 

The latter definition involves, besides the velocity field itself V = D°V , its gradient 
VV(r,£) = D 1 V. If the domain Q is bounded and the fluid fields belong to the NF set- 
ting both the scalar product and the norm can be trivially defined. However, this is not 
generally true in the case of an external (i.e., unbounded) domain Q, unless the velocity field 
V(r,i) vanishes for |r| — ► oo . This is usually denoted Sobolev- space functional setting (SB 
setting). 
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C. Open problems 



Referring to the 3D case the theory of existence and uniqueness of solutions is not yet 
complete in several aspects. 

In particular regarding existence, a major open question is whether strong solutions 
exist for all times (global strong solutions), i.e., the fluid fields A(r,t) remain at all times 
in the same (NF) functional class, or cease to exist in a proper sense, i.e., they develop 
spontaneously singularities which violate, at least locally, the Navier-Stokes equations. In 
both cases the answer would be of useful: in the first case to establish a global existence and 
uniqueness theorem with far-reaching consequences; in the second, to understand the nature 
of "singularities" or "irregularities" (i.e., for example, discontinuities arising in suitably 
higher-order derivatives of the fluid fields or violations of suitable bounds in appropriate 



functional spaces) possibly produced 



in a finite or infinite time interval 



Dy the spontaneous evolution and decay of turbulence 



23]. 



Another issue of fundamental importance in its own right is the problem of uniqueness of 
strong solutions, which are assumed to exist and to belong to a suitable functional class (for 
example NF or SB). Several authors have investigated the role of functional settings and 
supplementary regularity assumptions to be satisfied by the fluid fields and/or the data. In 
this regard, of fundamental importance is, in particular, the result obtained by Ladyzhen- 
skaya 19J, who was able to prove the uniqueness of the fluid fields A(r,t), assumed to exist as 
strong solutions belonging to the SB setting, by imposing the additional regularity constraint 
that the velocity field remains bounded in L 4 (Q). An analogous result was obtained later 
for strong solutions in L P (Q), with p > 2, (Kato and Ponce, 1984 26j; Sohr and von Wahl, 



1986 



Deuring and von Wahl, 1995 



301 ] ) and has been addressed by several authors 



also for weak solutions (see, for example, Iftimie, 1999 [32], Koch and Tataru, 2001 35], 
Galdi, 2002 [37]). However, these results fail, in general, in the case of an external domain 
Q, for the same reason indicated above. Therefore, a fundamental problem, as yet unan- 
swered, both for weak and strong solutions, is the search of the minimal functional setting 
for the existence of local and respectively global solutions in exterior domains (Galdi and 
Maremonti, 1986 j^]) and which do not decay at infinity in tt. Additional issues concern 
both uniqueness (J. Kato, 2003 [38J]) and existence (Giga et al.,1999 [33J; Giga et al. |34j). 
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D. Goal of the paper 



Goal of the paper is to address the problem of strong solvability of the incompressible 
Navier-Stokes equations (INSE) in Q x I. Q is here identified with a subset of the Euclidean 
space M 3 , formed at most by a finite number of disjoint, open and connected subsets of 
non- vanishing, and possibly non-finite measure (fluid subdomains); as a particular case, Q 
(fli) can also be identified with an unbounded subdomain of R 3 (external domain). Moreover, 
/ is subset of R, assumed either bounded or unbounded. In the following we intend to seek 
the minimal functional setting for the well-posedness of local strong solutions which, together 
with suitable derivatives of the fluid fields, in the case in which Q is an external domain, do 
not necessarily decay at infinity in Q (nondecaying strong solutions) and hence do not belong 
to Sobolev spaces. More precisely, this involves the search of appropriate minimal functional 
classes for the NS fluid fields V^r,^), p(r,t), assuring the validity of a local existence and 
uniqueness theorem for the initial-boundary value problem for NS equations. 

The treatment here developed is based on the IKT earlier developed by Tessarotto and 

nnn 

Ellero |3|, \% [5[, which relies on the construction of a suitable kinetic equation and providing 
exactly, as its moments equations, the required set of fluid equations. The IKT here adopted 
is shown to provide the basis for a local existence and uniqueness theorem for strong solutions 
of the NS equations. 

Key features of the present treatment follows: 

1) The fluid domain Q is a three-dimensional subset of M 3 which is not necessarily 
bounded, i.e., it can be an external domain. 

2) In the case in which Q is unbounded, the fluid fields, and in particular the fluid 
velocity V(r,t), are not required to be L P (Q) in the fluid domain Q, with p > 2. Hence, the 
fluid fields A(r, t) = {p,V,p}, together with their derivatives |i(r,t), Vi(r,i), V 2 i(r,t) 
and the volume force density f(r,t), are permitted to be nondecaying, i.e., to admit non- 
vanishing asymptotic values in Q. Thus, denoting e T the unit vector e r = r/r and B = 
{A(e r , t), §- t A, VA, V 2 A, f } , it must be such that 

lim rg n i | r |_ >OOJ B(r,t)= J B oo (e r ,t), 
0< t)| <oo, 

where the limits = {^(er, t), ^A^, VA^, V 2 ^^, f^} are allowed to be non-zero. 
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3) The fluid pressure p(r,t) is assumed non negative in the closure domain Q, i.e., p > 0, 
thus permitting the existence of vacuum regions, i.e., subdomains of M 3 in which p and p 
vanish identically, and corresponding free boundaries, i.e., the parts of (the border of) Q on 
which the fluid pressure vanishes identically. 

4) The boundary set dfl may generally be formed by piece-wise regular surfaces, curves, 
as well as isolated points of Q C R 3 . 

5) In validity of assumptions l)-3), a local existence and uniqueness theorem for strong 
solutions of the initial-boundary value problem for INSE is found by assuming that the 
fluid fields {p, V,p} and the volume force density acting on the fluid in Q x / satisfy 
stronger requirements than those provided by the NF setting and fl5ll6l) which are given by 
the following assumptions 

' {p,V,p,n,Q,f(r,t)} 6C(°)(0x/), 
V(r,t)GC( 3 > 2 )(ft x I), 
n,Q,p(r,t)GC( 2 ' 2 )(fi x J), 
f(r,t)eC( 2 -°'(nx I), 



(21) 



while 

[ V D (r)eC( 3 )(fi), 
n ) ,Q ,p (r)GC( 2 )(fi), (22) 

f(r,t )eC( 2 )(fi). 

Moreover, the velocity \~ w (r w (t),t) of each point, r w (t), of the boundary dfl is assumed at 
least of class 

V w {v w {t),t)eC^ 2 \l). (23) 

Here n(r,£), Q(r,£) and II o ( r ), Qo( r ) denote auxiliary fluid fields to be suitably defined. Fur- 
thermore, in the case in which VL is unbounded, the fluid fields A(r, t) — {p, V,p, n, Q} and 
the volume force density f (r,i), together with their derivatives J^(r, t), VA(r, t), V 2 A(r, t), 
are permitted to be nondecaying, i.e., to admit non- vanishing asymptotic values in Q. Thus, 
denoting e r the unit vector e r = r/r and B = {A(e r , t), j^A, V A, V 2 A, f } , it must be 



(24) 



where the limits B n 



lim re n i | r |_+ 00 .B(r,f)=B 00 (e r ,t), 
0< 1^(^,01 <oo, 

{A^ (e r , t), -j^Aao, VAqo, V^oo, f x: } are allowed to be non-zero. 



Assumptions (j2"Tl) - (l2"4"j) are denoted as minimal functional setting (MF setting). 
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E. Scheme of presentation 



The plan of the paper is as follows. First (in Sec. 2), the IKT developed by Tessarotto 



and Ellero 



flfl 



,|5| is recalled and its basic assumptions are pointed out. Basic consequences 
of the kinetic theory are analyzed in the subsequent sections 3 and 4. In particular, in Sec. 3 
the Navier-Stokes dynamical system, defined by the integral curves of the kinetic equation, 
is introduced and its conditions of existence, uniqueness and regularity are investigated. 
In Sec. 4, the problem of existence, uniqueness and regularity of the kinetic distribution 
function, for suitable initial and boundary conditions is analyzed, with particular reference 
to Maxwellian solutions. This permits us to obtain a theorem of existence and uniqueness 
of strong solutions of INSE in the MF setting. 

II. CONSTRUCTION OF IKT FOR INSE 

A remarkable aspect of fluid dynamics is related to the construction of inverse kinetic 
theories (IKT) for hydrodynamic equations in which the fluid fields are identified with suit- 
able moments of an appropriate kinetic probability distribution. Recently the topic has been 
the subject of theoretical investigations on the incompressible Navier-Stokes (NS) equations 
(INSE) Sflfl. The importance of the IKT-approach goes beyond the academic interest. 
In fact, fluid equations represent usually a mixture of hyperbolic and elliptic PDE's, which 
are extremely hard to study both analytically and numerically. As such, their investigation 
represents a challenge both for mathematical analysis and for computational fluid dynamics. 
For this reason in the past alternative approaches, based on asymptotic kinetic theories, have 
been devised which permit to advance in time the fluid fields, to be determined in terms 
of suitable moments of an appropriate kinetic distribution function. These methods, which 
approximate the fluid equations only in an asymptotic sense, are based on the introduction 
of suitably modified (fluid) equations which permit to advance in time the fluid fields only 
in an approximate sense. In particular, typically, their modified fluid equations actually 
describe weakly-compressible fluids. The discovery of IKT [1] provides, however, a new 
starting point for the theoretical and numerical investigation of hydrodynamic equations, 
since it does not require any modification of the exact fluid equations. In particular it holds 
for strong solutions, and permits to advance in time exactly the fluid fields by means of a 
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suitable kinetic distribution function /(x,i). Here x is the state vector x = (r, v), where 
respectively r and v denote the corresponding "configuration" and "velocity" vectors, and 
T is the phase-space spanned by x. In the following we shall assume that T is a phase-space 
of dimension 6. This is achieved introducing a phase-space classical dynamical system (here 
denoted NS dynamical system) 

x D -> x(t) = T t> t o x , (25) 

which uniquely advances in time the fluid fields by means of an appropriate evolution 
operator T tj t [3|. This is assumed to be generated by a suitably smooth vector field X(x, t), 

|x = X(x,t), (26) 
x(t ) = x G , (27) 

where x = (r, v) G T = fi x M 3 is a suitable state vector, r and v denoting suitable position 
and velocity vectors spanning respectively the configuration and velocity spaces f2 and R 3 , 
and finally the set of points x(t) (for t £ I) defines the (phase-space) Lagrangian trajectory 
of the NS dynamical system (NS Lagrangian trajectory). Therefore, introducing the corre- 
sponding pdf (probability distribution function) /(x,t) > 0, it fulfills necessarily in T the 
differential Liouville equation 

Lf(x,t) = Q, (28) 
where L is the Liouville streaming operator 

Eq. (1281) can be interpreted as a Lagrangian inverse kinetic equation. The corresponding 
equivalent Lagrangian form reads 

J(x(t), f)/(x(t), t) = /(x , t ) = / (x ), (30) 

where /(x(t),t) is the Lagrangian representation of the pdf, x(t) is the solution of the 
initial- value problem f[2"6"j) - f[2"T|) . / ( x o) is a suitably smooth initial pdf and 

dx(t) 



J(x(t),t) 

is the Jacobian of the map x D — > x(t). 



<9x a 



(31) 
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The vector field X(x, t) is in principle completely arbitrary. Therefore it can be defined 
in such a way that the inverse kinetic equation fl28l) satisfies an appropriate set of constraint 
equations and in particular so that it admits as a particular solution the local Maxwellian 
distribution 

/M(x,t) = ^3-exp{-^-}, (32) 

where u = v — V(r, t) is the relative velocity and v^r, t) denotes the thermal velocity 
defined in terms of the kinetic scalar pressure pi(r, t) (see below): 



2 Pl (r,t) 

v t h = \l ■ (33) 

Po 



As proven in Ref. 



thanks to the arbitrariness in the definition of the vector field /(x,i), 
and in the definition of the velocity moments of /, this permits us to construct an inverse 
kinetic theory for INSE. Hence, it follows in particular that: 

1. For prescribed initial pdf / (x ) the time-evolved pdf /(x(£), t), solution of the inverse 
kinetic equation (1281) . is uniquely determined by the NS dynamical system (|25|) . Thus, 
if the NS dynamical system, solution of the initial- value-problem fl26|) - fl27|) . exists then 
necessarily /(x(t),t) defined by Eq. (l30|) is a solution of the inverse kinetic equation 



2. For a prescribed choice of the vector field X(x, t), Eg. (1281) is an inverse kinetic equation 
for INSE, i.e., the fluid equations which define INSE - namely Eqs.([T])-(jl]) - are provided 
by a suitable set of velocity-moment equations of (12"5|) . Therefore, if /(x, t) is solution 
of the inverse kinetic equation (12"5|) . necessarily a suitable subset of its velocity moments 
coincide with the fluid fields {p, V,p} , which are solution of INSE. As a consequence, 
the corresponding velocity moments of Eq. (128]) must coincide identically with INSE. 



A. Assumptions of IKT 



The IKT approach for INSE and the corresponding NS dynamical system can be obtained 
in a straightforward way following the approach of Refs. ^, 3| • For this purpose, let us require 



that: 
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The pdf /(x,i) is summable in velocity space, in the sense that the velocity moments 



F G (r,t) 



d 3 vG(x,t)/(x,t) 



(34) 



exist in the closure domain Q x I and result suitably smooth in Q x I at least for the 
weight functions G(x,i) = l,v,E = |m 2 ,uu,£'u, where u = v — V(r,t) denotes the 
relative velocity (with respect to the fluid velocity). 



The fluid fields {p, V,p} coincide with the moments: 

p(r,t) = f v d 3 v/(x,t) = Po > 0, 
V(r,t)=±/ y d 3 vv/(x,t), 
p(r,t) =pi(r,t) - P , 
k PiM = f v d 3 vlu 2 f(x,t), 



(35) 



being P a a positive constant and Pi(v,t) denoting the kinetic pressure which is defined 
so that 



Pi(r,t) > P , (36) 

and hence it results p(r,t) > in Q x /. In addition we introduce the higher-order 
moments of the pdf: 

Q = [ d 3 vuEf, (37) 
Jv 

U = [ d 3 vuuf. (38) 
Jv 

Th 1S nnpUes that, by conduction, the ph ys ica, reahzabihty cond,tion s © and Q 
are identically satisfied. Furthermore, on can prove [3j that INSE are satisfied by 
suitably selecting the vector field X(x,i), i.e., by requiring that in the domain fl x I 
the moment equations 

I rf 3 vG(x,t) L/(x,t) = 0, (39) 
Jv 

corresponding to the first three moments, i.e., G(x.,t) = l,v,|u 2 , coincide with the 
equations of INSE. 

• Let us impose suitable kinetic initial conditions on the pdf /(x,t). The initial conditions 
are manifestly of the form 

/(x,* ) = / (x), (40) 
14 



where the initial pdf / (x) is assumed to satisfy the initial conditions (JTj),©,© for 
the fluid fields, which requires 

Po = Jyd 3 Vf (x), 

< V (r)=^/ v rf 3 vv/o(x), (41) 
k Po(r) = J v c/ 3 v|m 2 /o(x) - P , 

where u = v — V Q (r) and P > is an arbitrary initial constant. Instead, the initial 
moments 

Q = / d 3 vuEf (x), (42) 
Jv 

U a = I rf 3 fuu/„(x), (43) 

are arbitrary, so that - for example - it is always possible to impose that they vanish 
identically in Q [by suitably defining the initial kinetic distribution function / D (x)]. 



B. Additional assumptions - Extension of the NS Lagrangian trajectories on dfl 

An arbitrary Lagrangian trajectory defined by the NS dynamical system can gener- 
ally reach the boundary dQ, so that the solution of the initial value problem (j2"oT) - (l2"T|) 
x(£)= {r(£), v(£)} generally may not be defined in the whole existence domain of the NS 
fluid fields (J). This requires its extension on the boundary dQ. The result is achieved by 
defining suitable boundary conditions for the NS dynamical system, which must apply for 
an arbitrary phase-space trajectory x(t) of the same dynamical system. Certain restrictions 
must be placed on the possible motion of the boundary. This is due both to the assumed 
regularity of the fluid fields [i.e., the settings (I2T1) and (122)1 ] and the no-slip conditions to be 
imposed on the fluid velocity due to the Dirichlet boundary conditions (ITT1) . In particular, 
we shall require that: 

• If r w is an arbitrary point of dQ, its velocity, defined as V^(t) = V w (r w (t),t) = 
^rw(t), is by assumption a smooth real function of time in the sense (1231) . 

Nevertheless, the precise nature of dQ (i.e., if it is a surface, curve or an isolated point, 
remains in principle largely arbitrary. Thus, for example, dQ may be assumed as formed 
by piece-wise surfaces or curves, as well as isolated points, of M 3 . The definition of the 
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boundary conditions for the NS dynamical system can, in fact, be achieved in all such these 
cases. 

For this purpose, let us consider an arbitrary Lagrangian trajectory x(t) = {r(t),v(i)} 
which at time t = t c reaches the boundary dQ at the position rw(t c ) with nonvanishing 
relative velocity, i.e., is such that 

r(t c ) = r w (t c ), (44) 
hm \v(t) - V w (t)\ > 0. (45) 

Then, introducing the unit vector n w {t c ) defined so that 

«««c) = Km | r "'i;|" r '?' (46) 
t^t- \r w (t) -r(t)\ 

let us denote by x^(f c ) and x^(t c ), respectively, the incoming and outgoing Lagrangian 
trajectories, which are defined as: 



x (±) (£ c ) = lim x(t). (47) 



The boundary conditions for the NS dynamical system are obtained by imposing the bounce- 
back boundary conditions js| 

r(+)(i c ) = r(-)(t c ), (48) 
v (+) (t c ) - V w (t c ) = - [vH(t c ) - V„(t c )-] . (49) 



C. Kinetic boundary conditions 

To complete the set of assumptions required by IKT, appropriate kinetic boundary con- 
ditions must be defined for /(x,i). Consistent with Eqs. fl48l) and (|49p . they are achieved 
requiring that /(x, t) satisfies on dT the following constraints (A-C): 

• A) kinetic bounce-back condition: this is obtained by imposing the conservation of 
probability density at t c , i.e., 

/ (+ W(*c),v (+) (t c ),tc) = f { -\v w (t c ),2V w (t) - y ( -\t c ),t c ), (50) 

where ~Vw(t c ) = Vw( r w(tc),t c ) denotes again the velocity of the boundary dfl at the 
position r\y(t c ); 
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B) first fluid constraint: it is provided by the requirement that the first moment of 
the pdf yields the mass density, i.e., there results 



P 



[ d 3 vf(r w (t c ),v,t c ). (51) 
Jv 



C) second fluid constraint: it requires that, consistent with the Dirichlet boundary 
condition for the fluid pressure [provided by Eqs. ffTTj) ]. on dfl the kinetic pressure 
satisfies the constraint 

Pi(r w ,t) - P = p w (r w ,t), (52) 

where pw(jw,t) and Vw(rw,t c ) are prescribed in accordance with Eqs.l fTTi) and 
Pi(r w ,t), due to the position fl35|) . reads 



Pl(r w ,t c ) = J rf 3 v^{v- V w (r w (t c ),t c )} 2 f(r w (t c ),v,t c ). (53) 



D. Determination of the NS vector field 



Provided the fluid fields {p, V,p, n, Q} and the volume force density f(r, V, i) are con- 
tinuous in x / and suitably smooth, the vector field X(x,t) is found 0, [J to be of the 
form X(x,i) = {v, F} . Here the vector field F = F(x,£; /) can be written as 



F(x,t;/)=F (x,t;/) + F^x,*;/), 



(54) 



where the two vector fields Fq(x,£; /) and Fi(x,i;/), which depend functionally on the 
kinetic distribution /(x,i) (via the moments |p, V,p, II, Q}), are defined respectively as: 



F (x,t; /) = -V II - Vpi(r,t) + f + u • VV+z/V 2 V, 



(55) 



i r o 

Fi(x,t;/) = -u^—ln Pl - 



V 



Pi 



p-f-V+pV • VV - f - pV 2 V 

at 



+-v-Q--^[v-n]-Q 

Pi 2pi — 





f u 2 


-1} 









(56) 
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Then the following theorem (proven in Ref.[3() has the flavor of: 



THM. 1 - IKT for INSE 

Let us assume that: 

1) the fluid fields {p, V,p, n, Q} and the volume force f(r,£) belong to the functional class 
defined by the MF setting KM)-^): 

2) the pdf /(x,i) is strictly positive and is a particular solution of IKE [Eq. $28\)] in Qx I; 

3) the velocity moments of /(x,t) are defined by Eqs. (TffSP, and $3]) ; 

4) the pdf /(x,i) satisfies the initial conditions ( f^OP and 

5) the pdf /(x,£) satisfies on dVL the kinetic boundary conditions / f50j) . 137]) and / f5^j) . 
Tnen it follows that 

Tli) {p, V,p} are solutions of INSE in Q x I; 

TI2) i/ /(x, t) satisfies in the whole phase-space T = Qxl the initial conditions /(x, t Q ) = 
/m( x , t ), and on Q5 x I the boundary conditions 

f(+\r w (t c ),v,t) = f(-\r w ,2V w (t)-v,t), (57) 

£/ien it follows that /(x, i) zs solution of IKE [Eg. K28\) ] in the whole set T = Q x I. 

Tl 3 ) z/ /(x, i) is a particular solution of IKE, then {p, V,p} are solutions of INSE in 
ttx I. 

III. THE NAVIER-STOKES DYNAMICAL SYSTEM: A THEOREM OF LOCAL 
EXISTENCE AND UNIQUENESS 

Let us now analyze the existence, uniqueness and regularity of the IKT. In the following 
sections we intend to address in particular the well-posedness of the initial-value problem 
(!26|) - (l27|) which defines the Navier-Stokes dynamical system. 

It is immediate to prove that under suitable smoothness assumptions on 
{p, V,p, n, Q} , the problem (!26|) - (!27|) admits a local existence and uniqueness theorem 
and therefore defines a dynamical system Sns, to be denoted as Navier-Stokes dynamical 
system. 

THM. 2 - Local existence, uniqueness and regularity of the NS dynamical 
system 
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Let us require that the fluid fields {p, ~V,p,II, Q} and the volume force density f(r,t) are 
bounded in Q x I and belong to the MF setting K21\) -ft24\). 
Then it follows that: 

a) the vector field X(x,t) = [v, F(x,t)] is of class C < ^ 1 ' 1 \T x /); 

b) the solution of the problem j[2h)) -(2l\) exists locally and is unique: 

x (*) = X(* > *«» *) = M*)> v ( t )} = {Xr(* , t , t), Xv(x c , t a , t)} . 

Moreover: 

c) x(£) = x(x o ,t ,t) and its inverse, x Q = x(x(t),t,t G ) are at least: 

X(x ,t ,t),x(x(t),t,t ) G C^(T x / x /); (58) 

d) the Jacobian J(x(t),t) of the flow 

x -> x(£) = x(x , t , t) (59) 

is at least of class 

J(x(t),t) e c {1) (r x /) (60) 

and for all x Q G T and /or a// t a ,t G / zs non-vanishing and finite, i.e., 

J(x(t),t) ^ 0,oo. (61) 



Proof - In validity of the regularity assumption the vector field X(x,t) is manifestly of class 
(7'- 1 ' 1 )(r x J). Hence, thanks to the fundamental theorem of existence and uniqueness (see for 
example, Coddington and Levinson, 1955 Hirsch and Smale, 1974 [22 |) the regularity 
conditions (15"8~1) are implied. The additional result (1601 follows directly from Eq. ( 13TT) and the 
definition given above for the vector field F(x,£) [Eq. (l54l) ]. Moreover, thanks to Liouville 
theorem the Jacobian J(x(i),i) does not vanish or diverge for all x G G T and for all t , t & I. 
Q.E.D. 

An immediate consequence of THM.2 is the following 

COROLLARY 1 of THM.2 - Minimal functional setting of the fluid fields 

If the fluid fields obey the MF setting f21\) -(24\) there follows in particular that: 
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a) the NS vector field X(x,t) = [v,F(x,i)] results of class C^ 1 ' 00 ' 1 )^ x V x J); 
6j i/ie problem $2b)) -$ff\) admits one and only one solution x(t) 0/ c/ass (7( 1 > 2 ' 2 )(r x J x 7); 
cj i/ie Jacobian J(x(t),t) of the phase-flow is at least C^ 1 ' 2 ^ (F x J). 
Proof - The result follows from THM's 1 and 2 and the assumption of the MF setting 
for the fluid fields. Q.E.D. 

From THM.2 there follows the obvious further corollary: 

COROLLARY 2 of THM.2 - Extension of the solution x(t) on dtl 

The solution x(i) = x( x > ^o, of the initial value problem Eq. $MB) -$Wfy, prolonged on dQ 

by means of the bounce-back boundary condition f5W . exists and is unique. 

Proof - The result is an immediate consequence of THM.l, the assumption of validity of 

the MF setting and the boundary conditions previously introduced (see Sec. 2. 3) for the NS 

dynamical system. Q.E.D. 



IV. THE INITIAL-BOUNDARY VALUE PROBLEM FOR THE KINETIC PDF 
AND INSE: EXISTENCE, UNIQUENESS AND REGULARITY 

The previous theorem of local existence and uniqueness for the NS dynamical system, 
and its extension given by the 2nd Corollary of THM.2, can now be used to obtain a local 
theorem of existence, uniqueness and regularity for the kinetic distribution function /(x,£). 
Due to its arbitrariness, it is always possible to limit ourselves to the investigation of initial 
conditions of the form 

/(x,t ) = f M (x,t D ). (62) 

In view of the positions (l54l) . (|55l) (|56l) one can prove (see below) that this implies identi- 
cally /(x,£) = /jw(x,t). We stress that a more general result holds for an arbitrary (but 
suitably smooth and summable) distribution /(x,t) (see related discussion in Ref.jsJ). The 
choice (1621) warrants summability and existence of all the required velocity moments and in 
addition is consistent with the physical requirement set by the PEM (principle of entropy 
maximization; see Refs. 6j, l7|). Then, the following theorem holds: 

THM. 3 - Local existence and uniqueness of the kinetic pdf 

Let us require that: 
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1) the initial kinetic pdf / (x) coincides with the Maxwellian kinetic pdf /mOMo)/ 

2) the fluid fields {p,V,p}( r>t ) and the volume force density f(r,i) belong to the MF setting 
defined by Eqs.^)-(2^; 

3) the fluid fields satisfy the fluid initial and boundary conditions (0)-(GJ) and < f77]j ; 

^ /or a// (x,£) G T x / £/ie kinetic pdf /(x, t), £/ exists, satisfies the bounce-back 
condition and the conditions on pressure / f50j) . 137]) and 
It follows that in the domain T x I : 

a) /j^(x,t) zs a particular solution of IKE [Eq.fiStyl]: 

b) /jv/(x,t) differentiable; 

c) /jy(x,t) is summable; 

d) the velocity moments of ,/V(x,£), corresponding to G(x,t) = l,v,E = |n 2 , coincide 
with the fluid fields 

{p^,Pi=p + Po} {r , t) ; (63) 

and moreover: 

e) /]^(x,t) defined and summable also on the boundary set dT x J, where dT = dQ x V. 
Proof - The proof of a) follows by direct substitution of the position /(x,£) = /jvf(x,t) in 

Eq. (l28l) . In addition, Eq. (l30|) . shows that this solution necessarily corresponds only to the 
initial condition ( 1621) . i.e., there results necessarily in T x / : 

= 7f!|$ (64) 

b) Thanks to assumptions l)-4), /^(x,t) is manifestly differentiable in Tx I . c) Similarly, 
/m (x,i) is manifestly summable and d) its moments corresponding to G(x.,t) = 1, v,E = |w 2 
are by definition {p, *V,p% = p + -P }( r t y °0 Finally, thanks to Corollary 2 of THM.2, f M (x,t) 
is also manifestly defined on the boundary dVt. Q.E.D. 

THM.3 already contains in itself the basic ingredients required to reach the theorem of 
existence and uniqueness for INSE, which can be written as: 

THM. 4 - Local existence and uniqueness of the NS fluid fields 

In validity of THM.3 it follows that in the domain Q x J, the fluid fields {p, V,p}r T t) are 
necessarily strong solutions of the initial- boundary value problem of INSE. 

Proof - In fact, thanks to THM.l and 3 it follows that /m(x,£) is a particular solution 
of IKE [Eg. (1281) ] if and only if {p, V,^}^^ are necessarily strong solutions of the initial- 
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boundary value problem of INSE which belong to the MF setting. As a consequence the 
moments of fu (x,£), {p, V,pi}/ r ^ obey, by construction, INSE and the associated initial and 
boundary conditions, set by prescribing f M (x,t ) and respectively Eqs. (15Uj) . (15Tj) and (1521) . 
Q.E.D. 

We remark that, in the case Q is an unbounded domain, thanks to assumptions (EMI) , in 
principle no restriction is placed on the asymptotic behavior of the fluid fields for |r| — > oo 



m 
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25| . these solutions do 



at time t. Hence, in contrast to the customary approach 
not necessarily belong to Sobolev spaces. 

THM's. 3 and 4 contain the main contribution of the paper and the basic new results 
regarding the existence, uniqueness and regularity of nondecaying strong solutions of INSE. 
Thus, provided the fluid fields {p,V,p}r rt \ satisfy the assumptions of the MF setting in 
the set Q x I, together with the initial and boundary conditions, defined respectively by 
Eqs.flT])-© and (iTTj) . the solution of INSE exists and is unique. 

It is possible to show that the regularity assumptions on the fluid pressure are minimal 
in the context of the present approach. For example, if the fluid pressure is assumed to 
be unbounded from below, the inverse kinetic approach manifestly fails, since it cannot be 
related to the kinetic pressure in this case. In fact, the Maxwellian distribution is defined 
and results summable in velocity space (in the sense indicated above) only if the kinetic 
pressure is strictly positive pi(r, t) > 0. Moreover, if the partial time derivative ^ is assumed 
unbounded or discontinuous, the NS dynamical system cannot be defined any more, since 
the same NS vector field ceases to be at least continuous. 

The present result generalizes also the treatment of uniqueness of nondecaying strong 



solutions recently given by J. Kato 38j. The following remark is relevant in this context. 
First, there is a one-to-one correspondence between fluid fields and the Maxwellian distri- 
bution /m(x, t), while the uniqueness of /m(x, t) and of the phase-flow (143]) imply each 
other. Both are determined, in turn, by the solution of the initial- value problem (T26T) - (T2T|) 
and therefore by the value of the Jacobian J(x(£),£). Therefore, the uniqueness of strong 
solutions is intrinsically related to their existence. In fact, J(x(i),i) is defined and finite if 
and only if the fluid fields belong to the MF setting, and in particular the kinetic pressure 
Pi(r, t) results strictly positive, while - at the same time - the fluid fields are finite, i.e., 
Pi(r, £) < oo and |V(r, t)\ < oo. 
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V. DISCUSSION AND CONCLUSIONS 



In this paper the problem of existence of strong solutions of the initial-boundary value 



],y,3. 



problem of INSE has been addressed in the context of the IKT earlier developed 
The striking new features of present approach is that the fluid fields, and in particular the 
fluid velocity, do not necessarily decay at infinity in Q, hence the solution is not embedded 
in Sobolev spaces. 

The proof of existence, uniqueness and regularity of strong solutions of INSE is reached 
in two steps. The first one is based on the introduction of the NS vector fields X(x,t), 
suitably related to the fluid fields {p, V,p} , and by establishing the conditions of existence 
of its associated dynamical system, the NS dynamical system. The second step consists in 
the introduction of a suitable inverse kinetic equation and a related initial-boundary value 
problem, for which a theorem of existence, uniqueness and regularity can be reached. 

As a consequence, INSE are found to be expressed in terms of suitable velocity moments 
of the kinetic equation, whereas the initial and boundary conditions for the fluid fields are 
satisfied identically by proper definition of the initial and boundary conditions for the kinetic 
distribution function. 

A discussion on consequences and applications of IKT in fluid dynamics can be found in 
Refs.|g, 7| and 39 1 . 

Here we mention that by means of the IKT here adopted, the fluid pressure p(r,t) can be 
advanced in time self-consistently without solving explicitly the Poisson equation, since it 
can be determined directly as a moment of the kinetic distribution function by advancing in 
time the initial kinetic distribution function /(x, t Q ) in terms of NS dynamical system (|25|) 
and IKE ff28l) . As a fundamental result, an exact pressure evolution equation, advancing in 
time self-consistently the fluid pressure, can actually be achieved 401 ]. 

This feature is potentially relevant for the development of numerical solution methods 
for INSE based on the discretization of the kinetic distribution in phase space and yields a 
possible alternative to direct solution methods based on the solution of Poisson equation for 



the fluid pressure p 



3J (see also for example, 



41 



42. 



3,Q)- 
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